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It may prove useful in cosmology to understand the behav- 
ior of the energy distribution in a scalar field that interacts 
only with gravity and with itself by a pure quartic potential, 
because if such a field existed it would be gravitationally pro- 
duced, as a squeezed state, during inflation. It is known that 
the mean energy density in such a field after inflation varies 
with the expansion of the universe in the same way as radia- 
tion. I show that if the field initially is close to homogeneous, 
with energy density contrast \Sp/p\ <^ 1 and coherence length 
L, the energy density fluctuations behave like acoustic oscil- 
lations in an ideal relativistic fluid for a time on the order 
of L/\5p/p\. This ends with the appearance of features that 
resemble shock waves but interact in a close to elastic way 
that reversibly disturbs the energy distribution. 



I. INTRODUCTION 

The dark matter that in standard physics donrinates 
the mass of the universe might have originated in the 
same way as ordinary matter and radiation, in reactions 
at high redshift, or some or all of it might have been 
gravitationally produced by inflation ||l[, Q. A simple 
example of the latter is a single scalar field j/(x, t) with 
the pure quartic self-interaction potential 



V = A2/V4, 



(1) 



minimal coupling to gravity, and no other interactions. 
This could not be the dark matter that is gravitationally 
bound to galaxies and clusters of galaxies, but it might 
have an observationally interesting effect on structure for- 
mation. A nrodel based on the idea that the primeval 
fluctuations in this j/-fleld energy distribution broke ho- 
mogeneity to seed structure formation is compared to 
the measured power spectra of the radiation and matter 
distributions in l^ . Here I consider some aspects of the 
evolution of the j;-field energy distribution. 

Ford Q and Turner [g| show that after inflation ends 
the energy density in a homogeneous flcld with the quar- 
tic potential in Eq ([l]) scales with the cosmological ex- 
pansion parameter as py oc a(t)^^, the same as for ra- 
diation. This is important because it means the y-field 
energy could be large enough when inflation ends to be 
a signiflcant perturbation to the energy distribution, yet 
small enough thereafter that it does not interfere with the 
standard models for light element production at redshift 
z ~ 10^" and structure formation at lower redshift. It 
may be significant also that the parameter A in the quar- 
tic potential can be chosen so that eternal inflation gravi- 



tationally produces a potentially interesting value for the 
energy density py. There is the danger noted by Felder, 
Kofman and Linde that gravitational production can 
make py the dominant energy density, so this field as- 
sumes the role of the infiaton. Peebles and Vilenkin 0] 
show that in eternal inflation, with the quartic infiaton 
potential V^ = A0(?!)'*/4, this problem does not arise if the 
y-field parameter satisfies A^ ^ A. 

Since the mean value of py after inflation varies in the 
like radiation it nright not seem surprising that energy 
density fluctuations can behave in a similar way. An 
early step in this direction was taken by Khlopov, Mal- 
omed, and Zel'dovich M, who note that the quartic po- 
tential with positive A produces a positive effective pres- 
sure. Other steps toward the ffuid picture are reviewed 
in Section II. A linear acoustic wave model for the evolu- 
tion of small spatial fluctuations in the fleld energy den- 
sity can only be an approximation: the numerical solu- 
tion of Khlebnikov & Tkachev H shows the transfer of 
fluctuations in the spatial distribution of the energy to 
higher wavenumbers that Gruzinov B remarks can be in- 
terpreted as the relaxation to equipartition of the kinetic, 
gradient and potential energies. The purpose of this pa- 
per is to flnd a description of the onset of the departure 
from the acoustic wave model. 

This analysis assumes dy/dt = at the end of inflation, 
as follows from gravitational production of y, and the 
fleld value is close to homogeneous. Peebles and Vilenkin 
show this last assumption is a good approximation 
across the present Hubble length if A ^ 0.01. Khlopov, 
Malomed, and Zel'dovich |^ study the gravitational ef- 
fect of the evolution of the energy distribution in a scalar 
fleld with quartic plus quadratic self-interaction poten- 
tial, in a WKB approximation. Here I consider a pure 
quartic potential, where the effect of the peculiar gravi- 
tational acceleration is likely to be important only near 
black holes. I accordingly use an expanding spacetime 
that is homogeneous and isotropic. The y- fleld is sup- 
posed to have been squeezed by inflation to very large 
occupation numbers, so it can be treated as a classical 
fleld. 



II. PRELIMINARY REMARKS 

When the field is exactly homogeneous, the field equa- 
tion in Minkowski spacetime can be scaled to a one- 
dimensional oscillator, 



r] = -rj 



-1 < ?7 < 1, 



(2) 



where y = dy/dt, t is proper world time, and 'q{t) has 
unit aniphtudc. FoUowing Ford B, one arrives at a virial 
relation by multiplying Eq. (0) by rj and averaging over 
time. The left hand side is rji) — d{r]i))/dt — ff", in the 
average over time large compared to the oscillator period 
the second term is dominant, giving 0, 11^ 



(7)2) = (r;4) = 1/3, 



(3) 



where the energy is rf /2+ri'^ /A = 1/4 for unit amplitude. 
Greene et al. pO| show rj{t) has frequency 

uj = 2ti^I'^T{IIA)-'^ = 0M7 .... (4) 

The scaling to potential Ay*/4 and amplitude j/o is 

2; = J/o?/(Ai/2y„t), c^ = 0.847A1/V- (5) 

Turner Q notes that when the field is spatially homo- 
geneous, y = y{t), with amplitude yo, the time av- 
erage of the stress-energy tensor T/ is diagonal with 
p = y^/2 + A2/V4 = A(y4)/4 and p = {f /2 - Xy^/A) = 
A(yo/12), the time averages following from Eq. (||). The 
ratio pI p is the same as for an ideal relativistic fluid. 
One sees from the vanishing of the covariant divergence 
of this diagonal stress-energy tensor in comoving coordi- 
nates that the energy density scales with the cosmological 
expansion parameter as 



p ex a{ty 



(6) 



Other derivations are given by Ford B and Turner Iq] . 

If the field is inhomogeneous the action with quartic 
potential is simplified by changing from proper time t to 
conformal time t = J dt/a and letting y — ay. When lo 3> 
a/a this brings the action for y to the Minkowski form. 
In this limit a solution y{t, x) to the Minkowski problem, 
that has fixed mean energy density, produces a solution 
y{t, x) in which the kinetic, gradient and potential energy 
densities all scale as a{t)~'^ relative to the Minkowski 
solution. The mean energy density thus satisfies Eq. (0) 

(eg. i, 0). 

One gets some feeling for the evolution of the energy 
distribution in an inhomogeneous solution by starting 
from the homogeneous solution obtained by a Lorentz 
transformation of Eq. (J^) , 



y{r,t) ^ y^rj{X^^^y^-f{t 



r)). 



(7) 



The time average of the stress-energy tensor T^ in this 
solution is the same as for an ideal relativistic fluid with 
streaming velocity v, as one sees from the Lorentz trans- 
formation from Tij for y = y{t) or by computing Tij 
for y{r,t) in Eq. (m). If the field initially is locally well 
approximated by Eq. M), with amplitude yo and veloc- 
ity V that vary slowly with position, then Tij is close to 



the stress-energy tensor of an inhomogeneous relativis- 
tic ideal fluid, and the vanishing divergence of Tij tells 
us the energy starts to redistribute itself in the manner 
of a relativistic fluid. (This was pointed out to me by 
Vilenkin [12|). The analysis in the next section is meant 
to flnd indications of how long this fluid approximation 
might last. 

The discussion conveniently begins with yet another 
demonstration of the evolution of the mean energy den- 
sity in Eq. (0) when the field is inhomogeneous. The 
wave equation is 

d^y ^ ^ddy \7^y 



3- — 
a dt 



A2/^ 



(8) 



and the energy density and flux are 
p = 2/72 + (Vy) V2a2 + A2/V4, f = -yVy/a, (9) 

where the space derivatives are with respect to coordi- 
nate distance x and adx is a proper length. The energy 
conservation equation is 

dp/dt + V ■ f/a = -[3f + {Vyf/a^]h/a. (10) 

As in Eq. (m), the result of multiplying Eq. (H) by y 
and averaging over space and over a time interval large 
compared to the period uj~^ of oscillation of y and small 
compared to the expansion time a/ a is the virial relation 



(ij^) = {{Vyf)/a' + X{y'). 
This brings the mean energy density to 

(p)=3(2/2>/4+(Vy2)/4a2. 



(11) 



(12) 

The space average of Eq. (llOl) eliminates the energy flux 
density, and with Eq. (|l^) we get (p) ex a{t)~^. 

If the sum of the kinetic and gradient energy densi- 
ties is large compared to the potential energy then y has 
the familiar properties of a free massless field: initial 
large-scale fluctuations in the energy distribution grow 
smoother as in free streaming of a gas of relativistic par- 
ticles. In the situation to be discussed next the initially 
small gradient energy density grows larger in patches and 
may eventually become comparable to the kinetic and 
potential energies. In this case large-scale density fluc- 
tuations would evolve like a relativistic gas of massless 
particles with mean free path comparable to a character- 
istic de Broglie wavelength. 

The model presented in the next section is meant to 
describe the onset of this cascade to larger wavenumbers, 
when the fleld gradient energy density is small compared 
to the kinetic plus potential energy densities. The exam- 
ples in Section IV show that numerical solutions to this 
model are in satisfactory agreement with the energy dis- 
tributions from numerical solutions of the wave equation. 
Both demonstrate the failure of the linear acoustic wave 
model for the behavior of fluctuations in the energy dis- 
tribution, in the development of features that resemble 
shock waves except that they pass each other in a nearly 
free way. 



III. THE WEAK GRADIENT MODEL 



{{Vvf)/a^ - (3/4)/Vp. 



(17) 



Since the period of oscillation of y scales inversely with 
the amplitude (Eq. H), an initial gradient of the ampli- 
tude produces a growing gradient of the phase of the 
oscillation. This makes the field an oscillating function 
of position, as in Eq. (|^), and produces energy flux f that 
tends to carry energy from maxima toward minima of the 
distribution. When the field gradient is on the order of 
Vy ^ vy the streaming velocity in this rearrangement of 
the energy is ~ f / p ~ v. In acoustic waves in a relativis- 
tic fluid V is small when the density fluctuations are small, 
\5p/ p\ ^ 1. Following this condition, the weak gradient 
model (WGM) for the energy distribution in the y-field 
assumes two inequalities. First, the period ~ lo~^ of os- 
cillation of the field is much shorter than the wavelength 
~ {vu})~^ of spatial oscillations, because w <C 1. Second, 
the wavelength of field oscillations is much shorter than 
the coherence length L of the energy distribution. That 
is, the WGM assumes w » ww :^ 1/L. 

The goal is to study the redistribution of energy on the 
scale of L, not the more rapid oscillations of y in space 
and time, so the energy density and flux are smoothed 
over a window in spacetime that is large compared to 
the latter and small compared to L. The virial relation 
in Eq ( pT| ) applies to the averages over this intermediate 
scale, and the smoothed energy distribution accordingly 
satisfies Eq. (p^. The same smoothing of the energy 
Eq. (O) with Eq. (pi) yields 



(13) 



d{p)/dt + V • (f)/a = -4(p)d/a. 



The angular brackets mean the result of this smoothing. 
Eq. dl) is the first relation in the WGM. The second 
follows by writing out the time derivative of the energy 
flux density (Eq. H) and using the wave Eq. (g): 



dt 



4f- 



a V 4 



1 



;V 



1 



VyV^y. (14) 



Eq dl2) brings the mean of the first term on the right 
hand side to 

(Xy^/4 - yV2) = -{p)/3 {{Vyf)/{&a% (15) 

The last term is simplified by noting that the energy flux 
density satisfies 



'-->»v'»-^'^- 



(16) 



Eq. dig ) and the second term on the right hand side of 
Eq. ( ]15| ) will serve as small corrections to a wave equation 
for p and f , so I approximate these parts by their forms 
for a homogeneous energy distribution with streaming 
velocity v. Eq. (R) with v <^1 says Vy/a ^ — vy. With 
(y2^ ^ 4(p)/3 the energy flux density is (f) ^ 4(/9)v/3. 
Thus in the last term in Eq. (|l^) we have 



The last term on the right hand side of Eq. (nm is on the 
order of v'^ times the left hand side, and so is discarded. 
In these approximations Eq. (IIj) becomes 



ot a 



3a Sap Apa 



(18) 



Here and below I drop the brackets that indicate the 
smoothing of the energy density and flux. 

The conformal transformation that eliminated the ex- 
pansion of the universe in the action has the same ef- 
fect here. The transformations are p — a'^p, f = a'^f, 
di = dt/a. Without the tildes, Eqs. (^ and (|l|) in the 
new variables are 



dp 
di 

m 

dt 



-Vf, 



Ap 



-fV • f. 



(19) 



For motion in one space dimension this becomes 



dp 
'dt 



dx' 



Of 
dt 



IL. — __£. _jL — ___!. _ L—L 10()\ 



Idp 
3 dx 



f_dl 
p dx' 



with the wave equation 



dt^ 



1 d^f 



,2j-2 



I d^P 

3 dx^ 2p dx^ 



(21) 



Eq. ( |19D is the weak gradient model (WGM) for the 
evolution of the field energy density smoothed over scales 
small compared to the coherence length of the energy 
distribution and large compared to the wavelength and 
period of oscillation of the field. This is a computation 
in second order in the perturbation from a static homo- 
geneous energy distribution. 

The first part of the WGM in Eq. (|l|) just ex- 
presses energy conservation after the transformation to 
the scaled field in Minkowski spacetime. The second part 
can be compared to the behavior of an ideal fluid with 
rest energy density pr and pressure pr = Pr/i- The com- 
ponents of the stress-energy tensor T'^ for this fluid are 



00 



p - (472 - l)p,/3. 



r 



T°" ^ r = Apr-i^v°' /?,, 



(22) 



T 



Q/3 _ 



(5"/3 + 4^2^"^,/3)pj3. 



Here a labels the spatial Cartesian coordinates. The time 
part of the divergence of T*-' is the first part of Eq. (p^). 
The space part of the divergence, in second order pertur- 
bation theory, is 



at 



-Vp 



Ap 



v/^ 



-f(V.f)--(f.V)f. (23) 



If the flow is irrotational in first order perturbation the- 
ory, we can write (f • V)f = V/^/2 in the last term. 



In this case Eq. ( ]23| ) is the same as the second part of 
Eq.®. 

Since the evolution of the field energy distribution 
is the same as that of an irrotational relativistic fluid 
through second order in the perturbation from the static 
solution, the growing departure from the linear acoustic 
model may be expected to resemble the development of 
shock waves. One sees this in the WGM by considering 



The initial condition in Fig. 1 for a standing wave in 
one dimension is dy/dt = with 



solutions to Eq. (20) in second order perturbation theory. 



A standing wave solution is 

p — pi,[\ + A cos kx cos fit 
+ A^ cos2kx {\ ~ cos2Vtt - Vttsm2Vlt) /^ + . . ], (24) 

where pb is the mean energy density and Vt = fc/v3. A 
travelling wave solution is 



p ~ Pb[l + As.mk{x — f/\/3) 
- {ntA^/A) sin 2k{x - t/Vi) + 



(25) 



The nonlinear correction term in Eq. (E5^ makes the lead- 
ing part of a maximum of the energy distribution steeper 
than the trailing part, as in the development of a shock 
wave in an ideal fluid. In the standing wave solution the 
first nonlinear correction term vanishes when the leading 
term vanishes or is at a stationary point; at other times 
the correction term makes the extremum of the energy 
distribution more prominent and narrower. 



IV. NUMERICAL SOLUTIONS 

More detailed illustrations of the evolution of the en- 
ergy distribution, and tests of the WGM, are presented 
in Figures 1 to 5. The solid curves are computed from 
numerical integration of the wave Eq. (H) , and the dashed 
curves are from numerical solutions to the WGM in 
Eq. (poj). In these examples the initial field amplitude 
is almost constant at j/q. When the length and time unit 
is (A^'^2/o)~^, and the mean of y is scaled to unit ampli- 
tude, the field equation is 



d'^y/dt^ = V^y - y3 



(26) 



and the initial energy density is close to p = 1/4. 

The numerical integration of the WGM uses the den- 
sity as a function of position, rather than its mean value, 
in the denominator of Eq. dl3) , but that does not much 
matter because the energy density fluctuations remain 
small. The numerical solutions are periodic in the space 
intervals in the figures. The second spatial derivative is 
approximated as y'{ = (y^+i - 2yi + yi_i)/^^, where 5x 
is the spatial coordinate interval, with the usual exten- 
sion to two spatial dimensions. The time integration uses 
y{t + 5t) = y{t) + O.Ut[y{t + 6t) + y{t)], where y{t + 5t) 
is computed by iteration. 



y = 1 + B sm{kx). 



(27) 



The amplitude of the initial field fluctuation is i? = 
0.0125, so the maximum initial density contrast is Sp/ p = 
A~ AB ^ 0.05. The wavenumber is fc = O.OlA^/^j/o, or 
k = 0.01 in the units of Eq. (|2§). The labels in the fig- 
ure are the times from the initial condition in Eq. ( |27| ) 
in miits of the period T = 2-K\fijk of oscillation of the 
energy distribution in a linear acoustic wave model. The 
numerical solutions use 700 time steps per period of os- 
cillation of the field value at the mean amplitude yo-, or 
about 1 X 10^ time steps in the acoustic oscillation pe- 
riod T . Increasing the size of the numerical time step 
by a factor of two or four has little effect on the energy 
distribution at the time 20r shown in the figure. At 40T 
the numerical results are sensitive to the time step; ex- 
ploration of this stage of the evolution would require a 
more careful computation than has been attempted here. 
There are 350 space steps in the periodic boundary con- 
dition. If the coordinate interval is smaller than this, 
numerical instability spoils the integration before 20T. 

The wavelengths of the small-scale fluctuations behind 
the shock-like features in Fig. 1 vary with the coordinate 
interval 5x'i the value of the small-scale wavelength is a 
numerical artifact. I expect that in the limit of a contin- 
uous field the small-scale oscillations would approach the 
wavelength set by the period of oscillation of the field, 
since that and 8^ are the only relevant length scales in 
the computation. 

The top five curves in Fig. 1 show one half oscillation in 
the acoustic wave model. It is notable that the shock-like 
features pass through each other with little interference, 
leaving quite smooth energy distributions behind them. 

Fig. 2 shows the effect of lowering the initial ampli- 
tude B in Eq. ( |27| ) by a factor of ten and increasing the 
integration time by a factor of ten. As expected from 
Eq. (p4) , the time to appearance of significant departures 
from a linear acoustic oscillation scales inversely with the 
initial density contrast. It might be noted also that de- 
creasing B increases the wavelengths of the small-scale 
oscillations. 

Fig. 3 shows the energy distribution in a numerical 
solution to the wave equation in two spatial dimensions, 
for the purpose of checking whether the single spatial 
dimension in the first two examples significantly limits 
the nature of the departure from acoustic oscillations. 
The initial conditions are dy/dt — and 



y = 1 + B[svii{kxx) + sin(fcj,y)]/v2. 



(28) 



The last factor makes the rms density contrast the same 
as in Eq (|27|). The wavenumber on the horizonal axis in 
Fig. 3 is the same as in Fig. 1 , and the wavenumber on the 
vertical axis is 50/77 times the horizontal wavenumber. 



The energy contours are plotted at time 20T, where T 
is the acoustic oscillation period for the mode running 
in the horizontal direction. This is the same time as for 
the energy distribution at the top curve in Fig. 1, and it 
amounts to 12.99 times the acoustic period in the vertical 
direction. The amplitude B is the same as in Fig. 1, and 
the time step is the same, but there are only 50 x 77 
space steps. The energy distribution from the WGM is 
similar, though the density maximum at the lower left is 
narrower in the horizontal direction and the minimum at 
the upper right is broader. 

Better spatial resolution in this two-dimensional ex- 
ample would give a clearer picture of the failure of the 
acoustic wave model but would require a more ambitious 
computation. One does see the departure from acoustic 
wave oscillation, and it appears at about the same time 
as in the one-dimensional examples. There is no indi- 
cation in this example that the extra dimension allows 
qualitatively different behavior. 

In an application to a model for structure formation 
in cosmology , small-scale shock-like features the in y- 
field energy distribution could appear well before density 
fluctuations of interest to astronomy reach the Hubble 
length. Fig. 4 shows a check of the effect on the acoustic 
wave model for the large-scale energy distribution. The 
initial conditions are dy/dt = and 



y = 1 + _B[sin(fca;) — sin(8fca:)]. 



(29) 



The parameters in the numerical integration are the same 
as in Fig. 1. 

At time t = 2r, at two acoustic oscillation periods of 
the long wavelength mode in Fig. 4, the short wavelength 
mode has completed 16 oscillations and the cascade to 
shorter wavelengths is commencing. At i = bT the large- 
scale energy distribution is still close to the original sine 
wave, consistent with the acoustic wave model. This de- 
coupling from the nonlinear small-scale behavior is fa- 
miliar from other fluid models, and might be expected 
here on roughly similar grounds. In a region where the 
smoothed energy density is larger than average the y- 
field generally is oscillating more rapidly than average. 
That produces a gradient in the phase of the field os- 
cillation that produces a flux of energy away from this 
region. Local nonlinear features produce local fluctua- 
tions in the phase gradient but have little effect on its 
mean or on the mean energy flux density. 

Finally, a travelling wave solution offers a clearer il- 
lustration of the formation of shock-like features. The 
initial condition for the one-dimensional sinusoidal trav- 
elling wave in Fig. 5 is 



7;(x,i) = (l + 0)?7(t + (^), 
(j) ~ B cosk{x — t/v3). 



(30) 



to the wave Eq ( p^ if terms of order B^ and kB may be 
ignored. I use a numerical solution for r]{t), with rj = 1 
and r) = at x = in Eq. (^). This phase choice does 
not matter because rj oscillates much more rapidly than 
the energy distribution. The amplitude B and the other 
parameters in the computation are the same as in the 
standing wave example in Fig. 1. The growing departure 
from an acoustic oscillation looks very much like the de- 
velopment of a shock wave, as in Eq. (25). The time to 



development of large departures from the acoustic model 
is about half that of the standing wave, as might be ex- 
pected from the larger correction term in Eq. (p5|) and the 
coherence of phases of the leading and correction terms 
in the travelling wave. The length of the train of oscil- 
lations behind the shock-like feature doubles from time 
lOT to 20r. At 40T short wavelength oscillations in the 
energy density in the numerical solution fill space, in a 
complicated pattern of amplitude and wavelength fluctu- 
ations. This part of the evolution is not to be trusted, 
however, because the details are sensitive to the size of 
the space step in the computation. 



V. CONCLUDING REMARKS 

The discussion in Sec. II shows that under the initial 
conditions assumed here the initial departure from a ho- 
mogeneous energy distribution starts to rearrange itself 
in the manner of an ideal relativistic fluid with pressure 
Pr ~ Pr/i- Not so evident is that as the phase of the 
field oscillation as a function of position winds and un- 
winds the energy distribution continues to behave like a 
fluid. This behavior follows in the weak gradient model 
in Sec. Ill (Eq. 11^). The weak gradient model in turn is 
equivalent to the equation of motion of a relativistic fluid 
to second order in the perturbation from a homogeneous 
energy distribution, provided the flow is irrotational in 
first order perturbation theory. 

In the numerical tests in Figs. 1 to 5 the weak gradient 
model agrees well with the energy distribution obtained 
from numerical solutions of the wave equation. That is, 
we have good evidence for the picture of how the near 
acoustic wave behavior of the energy distribution ends: 
features that resemble shock waves appear after a time 
on the order of 



r„i - T/5p/p, 



(31) 



where rj{t) is a solution for the quartic oscillator with 
unit amplitude in Eqs. (|^) and (||). Eq. (|3^) is a solution 



where T is the oscillation time in the acoustic model and 
the rms density contrast is 5p/ p. 

The features are not true shock waves, of course. In 
particular, they excite energy fluctuations in a reversible 
way: the field returns to a smooth energy distribution 
after the feature has passed by. There is no artificial 
viscosity in the computation, and I have not been able 
to see how the effect could be a numerical artifact; this 



seems to be real reversible process. The near elastic in- 
teractions illustrated in the half cycle in Fig. 1 remind 
one of solitions, but this behavior has a limited lifetime: 
the length of the train of oscillations behind the leading 
edge increases with time. 

More ambitious analyses than have been attempted 
here would be needed to explore the evolution of the 
shock-like features. If, as seems likely, the fluctuations in 
the spatial distribution of the field energy eventually cas- 
cade to noise extending to wavenumbers comparable to 
the frequency uj of the field oscillation (Eq. ^ , the energy 
distribution may be expected to end up behaving like a 
gas of interacting particles with de Broglie wavenumber 
on the order of w [^ . The near fluid behavior on scales 
much larger than the de Broglie wavelength might still 
obtain, however, because the gas has a short mean free 
path. This aspect of the evolution remains to be ana- 
lyzed. 
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FIG. 1. Evolution of a one-dimensional standing wave in 
the energy distribution of a field with a quartic self-interaction 
potential. The solid lines are from a numerical integration 
of the wave equation (Eq. Kq), and the dashed lines are a 
numerical solution of the weak gradient model (Eq. EG). The 
energy distributions are plotted at the indicated multiples of 
the period of oscillation of the initial sine wave in the linear 
acoustic wave model. 
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FIG. 2. Illustration of the effect of decreasing the wave 
amplitude in Fig. 1 by a factor of ten. The curves are plotted 
at the indicated multiples of the acoustic oscillation time. All 
other parameters are the same as in Fig. 1. 




FIG. 3. The two-dimensional energy distribution after 20 
times the period for an acoustic osciUation in the horizontal 
direction. This is 13 times the acoustic oscillation period 
in the vertical direction. The energy distribution is from a 
numerical integration of the field equation with the initial 
condition in Eq. (Eq). The time step, amplitude parameter 
B, and wavenumber in the horizontal direction are the same 
as in Fig. 1, but the spatial resolution is worse by a factor of 
seven. The highest density contours, in the two peaks at the 
lower left, are at contrast 5p/p = 0.069; the lowest contour, at 
Sp/p — —0.045, is at the upper right. The nearly flat regions 
in the upper left and lower right are at close to zero density 
contrast. 
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FIG. 4. Illustration of the effect of small-scale nonlinear 
fluctuations on evolution of the mass distribution on larger 
scales. The initial condition is in Eq. (E9|). At the time of the 
top curve the long wavelength component has oscillated five 
times and the short wavelength component 40 times in the 
acoustic wave model. The other parameters are the same as 
in Fig. 1. 
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FIG. 5. A travelling wave solution with the initial condition 
in Eq. (|3Cl). All other parameters are the same as in Fig. 1. 



